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Abstract. Assuming the form of the entropic dark energy as arises form the surface 
term in the Einstein-Hilbert’s action, it’s evolution were analyzed in an expanding flat 
universe. The model parameters were evaluated by constraining model using the Union 
data on Type la supernovae. We found that in the non-interacting case, the model 
predicts an early decelerated phase and a later accelerated phase at the background 
level. The evolution of the Hubble parameter, dark energy density, equation of state 
parameter and deceleration parameter were obtained. The model is hardly seems to be 
supporting the linear perturbation growth for the structure formation. We also found 
that the entropic dark energy shows phantom nature for redshifts z < 0.257. During 
the phantom epoch, the model predicts big-rip effect at which both the scale factor of 
expansion and the dark energy density become infinitely large and the big rip time is 
found to be around 36 Giga Years from now. 
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1. Introduction 

The mounting observational evidences [H [2] have proved that the current universe is 
accelerating in expansion and the acceleration was began in the recent past in the 
history of the universe. It was conjectured that this acceleration is caused by the 
exotic form of energy called dark energy (DE) which can able to produce negative 
pressure, unlike the conventional form of matter. The origin and evolution of dark 
energy is still the most important problem for the physicist today. The standard ACDM 
model in which the universe consists of cod dark matter and cosmological constant as 
dark energy is matching well with the observational data, but it is facing the severe 
cosmological constant problem that the theoretical value of the dark energy density 
as predicted form the held theoretical considerations is nearly 120 orders larger than 
the observational value [3111]. This motivates the consideration of the dynamical DE 
models. Various DE models were proposed aiming mainly on solving the cosmological 
constant problem. For a review of dark energy models, see, Copeland et. ah [5], Sahni 
et. al. [3] and references therein. Among the different classes of models the quintessence 
models [SIEI and K-essence [H] models based on the scalar helds, while Chaplygin gas 
model[9] is an example for the one trying unify both dark matter and DE. Various 
other models are existing in the current literature, for instance, holographic dark energy 
models [ini, Agegraphic[II] models etc. Among these variety of models, entropic-force 
DE model gained much attention recently. Entropic DE (EDE) model was hrst proposed 
by Easson et. al. da 113]. This model is based on the idea of entropic gravity, proposed 
by Verlinde|14j. Verlinde explained gravity as a kind of force related to the change in 
entropy, hence the name entropic gravity, and he derived the held equations of gravity 
from the second law of thermodynamics. The entropic dark energy in the entropic 
gravity model is being arised from the surface term, often neglected, appearing in the 
Einstein-Hilbert’s action. It was conjectured that the horizon, the boundary of the 
universe is acting as the so called ‘holographic screen’. 

In the work of Easson et.ah, it was shown that surface part in the action would 
add a positive term to the acceleration equivalent to ChH^ + CgjH, where H is the 
Hubble parameter, Ch^C^ are the model parameters assumed to be in the range 
3/27r < Ch < 1 and 0 < < 3/27r, where over-dot represents a derivative with 

time. They have demonstrated from error plot analysis using Type la supernovae data, 
that, the model is in good agreement in predicting the distance modulus of various 
supernovae and the universe moves smoothly from a decelerating to an accelerating 
epoch at around a redshift z ~ 0.5. Later this model was considered by Spyros et. 
al.[T5], as a running vacuum, with constant equation of state, -1 so that the model 
leads to either eternal acceleration or eternal deceleration for the respective parametric 
range hence doesn’t predict a transition from an early decelerated phase to a latter 
accelerated phase of expansion and hence they argued that EDE model is not viable 
both at the background level and perturbation level. For obtainig this result they have 
assumed a general conservation law, which implies a energy transffer between the EDE 
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and other cosmic components. On the other hand Easson et ah doesnt assnme any 
energy transffer, hence each component assnmes sepertate conservation law. 

In the present work we have assnmed the general form for the EDE density and 
derived the Hnbble parameter by following the Easson et al. approach, where there is 
no interaction between the dark energy and other cosmic componens, especially matter. 
The model is then contrasted with the observational data on Type la snpernovae to 
extract the model parameters. The various cosmological parameter like dark energy 
density, equation of state, deceleration parameter were calculated and their evolution 
are studied. The paper is organized as follows. In section two we obtained the Hubble 
parameter and constraint the model with observational data to extract the model 
parameters. The error bar plots with snpernovae data were constructed and with the 
model prediction. The model is also compared with the Stern et al. data on Hubble 
parameter at various redshifts. Section three consists of our analysis of the various 
cosmological parameters and also results on the Om diagnostic analysis of the model. 
In section four we have presented our observations on the big rip singularity occurred 
in this model followed by conclusions in section Eve. 


2. Entropic dark energy, Hubble parameter and estimation of model 
parameters 


The surface term in the action, which is often neglected, will causes the EDE. The 
Einstein Hilbert action including the surface term can be schematically expressed 

as[T^ ITB] 

1= [ {R + Lrn) + ^ I K (1) 

JM JdM 

where R is the scalar curvature, is the Lagrangian corresponds to matter, hied and K 
is the trace of the extrinsic curvature of the boundary. The variation of this action with 
respective boundary contribution will lead to the Einstein’s held equation, which when 
combined with Friedmann metric, gives the Firedmann equations, of which especially 
the acceleration equation become |T2], 

^ = + + ( 2 ) 

In this equation the last two positive terms on the right hand side constitute the entropic 
dark energy density, with coefficients and Cr, which constitute the EDE density. 


Pede — 3 ( CrH + Cjj 


H 


(3) 


where Cr and are the model parameters and we follow the natural system of 
universe, SirG = c = 1. For a spatially hat Friedmann universe, which is being favored 
by cosmic microwave background (CMB) observations and predicted by inhationary 
models [niiiH], this DE satishes the Friedmann equation. 


j^2 _ pm I PEDE 

~ ” 3 ” 3 


( 4 ) 


3 
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where pm is the matter or effectively cold dark matter density, and we nse the natnral 
system nnit in which SttG = c = 1. Assnming no interaction between dark energy and 
non-relativistic matter, the matter and entropic DE satisfies the conservation eqnation 
separately, 


Pm + {pm + Pm) — 0 

(6) 

Pede + {pEDE + Pede) = 0 

(6) 


where Pm/EDE are the pressnres of the matter/EDE components of the nniverse. 

Eqnations (|^ and Q were properly combined to obtain the differential eqnation 
for Hnbble parameter evolntion, 

+ + = o ( 7 ) 

dx \ Cp J Cp 

where h = H/Hq, Qmo is the present valne of mass parameter of dark matter and 
X = logo with a as the scale factor of expansion. The solntion of this will gives the 
evolntion of Hnbble parameter as. 


h‘^{x)=rie {1 - p) e b/CjiA 


( 8 ) 


where. 


V 


n 


mO 


l + ilc 


H 



(9) 


which shows that p = flmo if Ch = ^Cp/2, p < Qmo if C'h < 3C^/2 and p > Hmo 
if Ch > 2>Cp/2. The Hnbble parameter given by the above eqnation also satisfies the 
condition h = 1 for z = 0 corresponds to x = 0, apt for the cnrrent state of the nniverse. 
Especially for the case Ch = ‘iCp/2 the Hnbble parameter rednces to, 

+ (1 - , ( 10 ) 


and with additional constraint Hmo = 1) corresponds to the matter dominated case the 
Hnbble parameter become, 

( 11 ) 


which in fact the Einstein-de Sitter model. 


2.1. Contrasting with Type la data 

We have evalnated the model parameters Ch, Cp and also the present valne of Hnbble 
parameter, Hq by constraining model nsing the observation data. These parameters are 
determined by minimizing the fnnction, 

^2 _ ^ [Tobsi^i) ( 12 ) 

i 

where Pobs{,Zi) is the observed distance modnlns of the snpernova at redshift 2 ;*, pth{zi) is 
the corresponding theoretical modnlns and Uj is the nncertainty in the SNe observation. 
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Model 

A. min 

xLn/d-O.f 

a 


Ho 

EDE with klmo = 0.30 

311.29 

1.024 

0.813 

0.415 

70.59 

EDE with Dmo = 0.27 

311.29 

1.024 

0.832 

0.374 

70.59 

EDE with Dmo = 0.25 

311.29 

1.024 

0.844 

0.346 

70.59 

ACDM model 

311.93 

1.026 

- 

- 

70.03 


Table 1. Best estimates of the parameters using supernovae data for different 17^0- 
The term, d.o.f = degrees of freedom = N - n, N=307, the number of data points, 
n=3, number of parameters in the model. For comparison we have also evaluated the 
corresponding values for the standard ACDM model for the same data set. 


We have used the “Union” SNe data set [19], which consists 307 type la supernova. The 
theoretical distance modulus is, 


Hth{zi,a,(3,Ho) = Slog 


10 


dj 


Mp(^ 


+ 25, 


where is the luminosity of the supernova and is given by, 

dz' 

dL{zi,a,(d,Ho) = c(l + z) 




(13) 


(14) 


where c is the light speed. The best estimates of the parameters Ch^C^ and Hq for three 
different values of flmo = 0.25, 0.27,0.30 are obtained by minimizing the function and 
are shown in table The x^j„=311.29 and the Hubble parameter Hq = 70.59 are the 
same for the three values of the dark mass parameter but Ch parameter value decreases 
as Hmo increases and parameter value increases with Hmo- For a comparison we 
also evaluated the Xmin for the ACDM model, which give values almost close to the 
EDE values. The statistical correction for the parameters values were obtained by 
constructing the conhdence regions in £gure[^for 11^0 = 0.3. The best £t corresponds 
to 99.73% are Ch = 0.813 ± 0.056, = 0.415 ± 0.061 and for 99.99% probability are 

Ch = 0.813 ± 0.072, C^ = 0.415 ± 0.079. 

The best estimated values shows that 2{Ch — 1)/C^ ~ —0.902, which implies that 
at large redshift, z >> 1, the Hubble parameter become, 

~ r]{l + zY, (15) 


ensures that the early universe has gone through a matter dominated phase. In 
comparison with the corresponding equation in standard ACDM model, i.e. = 

f^mo(l + z)^, the parameter rj is now corresponds to the dark mass parameter Dmo 
in ACDM model, however it is slightly large, rj = 1.23flrno, as per the best estimates 
of the model parameters. This hike is due to the fact that the EDE is mimicking the 
dark matter characteristic in the early matter dominated epoch, a fact to be clear from 
the considerations of equation of state parameter, and is contributing towards the total 
non-relativistic mass density at z >> 1. Similar situations of increased mass parameter 
were discussed in many recent literature, for instance the mass parameter = 1 . 111^0 
in the quintessence model considered in reference [22] . However for the sufhcient growth 
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Figure 1. Confidence intervals for the parameters a = Ch and (3 = for the 
Hubble parameter Hq = H{z = 0) = 70.42 and tlrno = 0-3 using the Union supernovae 
data. The confidence intervals shown are corresponds to 99.73% probability (inner 
one) and 99.99% probability (outer one). The dot represents the values of {Ch,C^) 
corresponds to xLin- 



Redshift 


Figure 2. Comparison of the distance modulus between the EDE model with 
= 0.3 and supernovae data. The error bars are corresponds to the observational 
data and the continuous line is the prediction from the present model. 


of perturbations in the matter dominated epoch of a universe consists of quintessence 
filed, it is argued that flmo ^ 1-15^2^0 |23]. Taking this constraint in to consideration, it 
seems that the non-interacting EDE model hardly supporting the perturbation growth, 
required for the structure formation in the matter dominated era, hence it may not be 
viable at the perturbation level. 

Using the best estimates of the model parameters, we have compared the EDE’s 
prediction of the distance modulus di, of supernovae at various redshift in equation (14) 
with supernovae data and is shown in figureThe figure shows good agreement between 
theory and observation at low redshifts. Similarly the Hubble parameter evolution in 
EDE model is compared with the Stern data [21] and also with the ACDM model in 
figure It is seen that the Hubble parameter evolution of the EDE model is htting 
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Red shift 


Figure 3. Hubble parameter evolution of EDE model with Vtmo compared with Stern 
data on Hubble parameter and ACDM model, the error bar lines are corresponds to 
Stern data, dotted line corresponds to Hublle parameter evolution of EDE model and 
continuous line corresponds to ACDM model. 

fairly well with the observational data at lower redshift but high redshift the htting is 
not seems to be good. At the low red shift region for the expansion rate in the EDE 
model is slightly lagging behind the the ACDM model, at around z ~ 0.7 the expansion 
rate in the non-interacting EDE model making cross-over into large rate and is very well 
leading ACDM model at large redshift. This is due to the fact that at very large redshift 
redshift the EDE dark energy is almost having the behavior of cold dark matter. While 
discussing the model independent extraction of dark energy properties, Sahni et. ah 
PU] shows a similar type of contrast in the behavior of the Hubble parameter between 
the model independent dark energy and ACDM model. 

3. Evolution of cosmological parameters 

In the present model, the dark energy density evolution with respect to redshift can be 
obtained form equation (|^ as, 

klEDE = he"'" + (1 - h) (16) 

In terms of the scale factor the above equation become, 

klEDE = {v ~ -b (1 — h) _ (17) 

For the best estimates of the model parameters, it can be easily seen that 2{Ch — 
1)/C^ — so the last term in the above equation become proportional to a^. Hence 
in the remote past of the universe (at large redshift, 2 ; >> 1) corresponds to a —)■ 0 the 
EDE density become, ELede — iv ~ which corresponds to the behavior similar 

to the cold dark matter. But at small redshift corresponds to large a the dark energy 
density, Qede {1 — rj) a, grows almost in proportion to the scale factor. 

The equation of state parameter can be obtained from standard relation |25] in 
terms of redshift as, 

1 /3(h - D™o)(l + ^)' + - h)(l + \ 

~^3 {r] - Qrno){l + zf + {1 - r]){l + j 


( 18 ) 
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redshifl 


Figure 4. Evolution of the equation of state parameter for the best ht values of the 
parameters. 


At very large redshift, z » 1 only the terms proportional to (1 + will effectively 
contribute, and other terms will be negligibly small, as a result ojede —^ 0. That confirms 
that in very early stage the EDE has a behaviour similar to the cold dark matter. But 
as 2 ; —>■ — 1 the equation of state will behaves as, lvede —)■ —1 + 2{Ch — which is 

for the best estimates of parameters (for Q^no = 0.3), become, u)ede{z —>• —1) — >■ —1.3. 
The evolution of the equation of state parameters given hgure The figure shows that 
EDE equation of state parameter crosses the phantom divide, ojede = —1, and proceed 
towards still lower values in the latter evolutions, before stabilizing at around —1.3. 
During it’s phantom behavior at which the equation of state is less that -1, the EDE 
violates the null energy condition, {p + p) >0 and the energy density increase during 
the further evolution. In fact we have already seen form the previous paragraph that 
the EDE density does shows such a behavior in the late universe. From literature one 
can End that phantom behavior can be constructed from scalar field with a negative 
kinetic term, 

Pde =+ V{(j)) Pde = (19) 

which will leads to 


UlDE — 


-u^ + v{<fy 


( 20 ) 


Where 0 is the scalar held, V (0) is the potential energy and the above equation gives 
(^DE < — 1 if 02 > 0. The current value of the equation of state in the present model 
is noted as ujede{.z = 0) ~ —1.185. This values is higher than the observationally 
constrained value, around ojede{z = 0) ~ —0.93, |2H] obtained form the joint analysis 
of WMAP+BAO+SNe data. The present value of the equation of state parameter of 
this model indicating the phantom nature of the entropic dark energy in the current 
phase of the universe. It is to be noted that some latest results on the equation 
of state in fact prefer phantom dark energy models. For example one of the latest 
cosmological data gives, ojde = l-04j^Q;5Q[26], while in reference [2H], the equation of 
state parameter is deduced as ojee = The recent results from Plank satellite 

gives, ude = -1.49l°|^ 
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Redshift 


Figure 5. The evolution of the deceleration parameter with redshift. 


A forecast of the phantom behavior of dark energy is the occurrence of big-rip 
effect or future singularity [27115I]| 151] . due to the tremendous increase in the dark energy 
density as the universe expands. We have already shown that the density is increasing 
in proportion to the scale factor at sufficiently lower redshifts. So the big rip is seems 
to be inevitable in this model. We will analyze this fact in a later section. 

The deceleration parameter can be obtained using the relation, 

1 dh"^ 



Using the expression for from equation([^, the evolution of the deceleration parameter 
is, 


9 = - 1+2 


3))(1 + zf + - )))(! + 

r/(l + z)^ + (1 — '?)(! + 


( 22 ) 


In the limit z » 1 the second terms both on numerator and denominator inside the 
parenthesis on the right hand side of the above equation will be vanishingly small and be 
neglected, as a result, g —)■ 1/2, con&ming that the EDE behave like CDM in the earlier 
period[]. More over at the special condition, 3C^/2 = Ch, = 1 at which rj = 1, the 
deceleration parameter become q = 1/2 implies Einstein-de Sitter universe as we already 
noted in the previous section while dealing with the evolution of the Hubble parameter. 
As the universe evolves further the deceleration parameter deceases and become negative 
at the recent past in the evolution of the universe. The negative value of the deceleration 
parameter implies acceleration in the expansion of the universe. For small redshift, 
when the EDE would increasingly dominate, the deceleration parameter be stabilized 
at g —)■ — 1 -|- {2Ch — 1)/C^, which is in fact less than —1 as per the best £t values of 
the parameters Ch and C^. This behavior of the parameter indicates a transition from 
an early decelerated epoch to a later accelerated phase of expansion, and this transition 
is found to be occurred at a redshift of zt = 0.57. Figure shows the behavior of the 
deceleration parameter with redshift. The present value of the decelartion parameter 
is, q{z = 0) = —0.8. The observational constraints on the parameters, g and zt are 
q{z = 0) = —0.64 ± 0.03 [32l l33]and zt = 0.45 — 0.73[3l]. Even though the predicted 
value of transition redshift zt in this model is almost in agreement with the observational 
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0.0 0.5 1.0 1.5 2.0 

scale factor 


Figure 6. Evolution of DE with scale factor. The plot shows that when the DE 
density is decreasing in the early stage but increasing at the later stage where the 
equation of state is less than -1. 


constraint, the present value of the deceleration parameter is not in exact agreement 
but slightly less than the observational value. 


4. Big rip singularity 


In the earlier section we have seen that the EDE shows phantom nature for 2 ; < 0.257. 
During phantom phase dark energy density ELede will increases in proportion to the 
scale factor, see hgure In this section we want to check whether the phantom nature 
of EDE indicates a big rip scenario. In the phantom stage the universe is increasingly 
DE dominated, hence the Hubble parameter given in equation([^ can be approximated 
with the second term on the right side alone, 

~ (1 - T]) exp(-2(C// - 1)x/Ch). (23) 


Considering the best £t values of the parameters, Ch and it can shown from the 
above equation that, the scale faactor will diverge in hnite duration of time in future, 
gives as big-rip time, Uip, 


f ■ 

^rip — 


(1 


(24) 


From the equation ([^ of rj, the big rip time t^p is strongly depending on [2Z]- 
With Vtmo = 0.30, the above equation leads to trip — ^0 — 2.5Hq^ for the best estimates 
of the model parameters, and is around 36 Giga Years for present model. While EDE 
density behaves as klede oc a in DE dominated phase, it also will diverge during the big 
rip. The big rip phenomenon due to the phantom nature of the dark energy was first 
discussed by Caldwell|27j. The dark energy effectively anti gravitating in nature. In 
the phantom phase, the effect of this anti gravity is so high that it dissociates the large 
scale structures [27] in the universe. During this stage the the expansion rate H grows 
with time, see figure]^ so that the horizon (whose size is proportional to H~^) closes in 
on as and the galaxies will thus crosses the horizon to disappear. 
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redshift, z 


Figure 7. Growth of the expansion rate especially for z < 0. Compared to ACDM 
model (the dotted line) the expansion is growing at a faster rate when z < 0. 


5. Conclusions 

Entropic dark energy (EDE) was proposed first by Easson et al.[T2] based on the 
Entropic gravity model. Technically the EDE arises from the surface term in Einstein- 
Hilbert action of gravity, and having a general form ChH^ + C^H. Easson et ah have 
shown that the model predicts a transition form a decelerating phase to accelerating 
one at around a redshift, = 0.5. They have also found a good agreement between the 
theoretical distance moduli of supernovae predicted by the model and the corresponding 
observational moduli. We have analyzed the non-interacting EDE model, in which both 
the dark energy and dark matter satishes seperate conservations laws, by evaluating the 
best estimates of the model parameters using the Union data on Type la supernovae. 
After deriving the Hubble parameter, the evolutions of the dark energy density, equation 
of state parameter and deceleration parameter were studied. The model predicts an 
early decelerated epoch and a later accelerated epoch. The model will reduces to 
the Einstein-de Sitter universe for the special conditions of parameters, 3C'^/2 = Ch 
and Dmo — 1- The model predicts a slightly higher value for the mass parameter 
around, Qmo — l-25Dmo, which may dilute the linear perturbation hence may negatively 
affect the structure formation. However a hrm conclusion regarding this needs detailed 
calculation. The evolution of the deceleration parameter is studied and found that there 
is a transition form the early decelerated epoch to a later accelerated epoch at around 
a redshift z = 0.57, a value slightly higher than that predicted in the work of Easson at 
al., but well within the WMAP range [25]. 

The equation of state is found crossing the phantom divide at ^ < 0.257 and 
stabilizes at around —1.3. It’s present value is around uede ~ —1.185. We have also 
shown that the phantom nature of the EDE leads to big rip situation, such that during 
a hnite time in the evolution of the universe the scale factor and dark energy density 
were become inhnitely large. The big rip time is found to be around 36 Giga years from 
now. As pointed out by Caldwell, this may leads to the dissociation of structures and 
apart from that the horizon will approach us hence the structures may cross the horizon 
to disappear. 
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However one may the note the work of Spyros et ah na at this juncture, where the 
authors treated the entropic dark energy as a dynamical vacuum with constant equation 
of state equal to -1, so that the dark energy is interacting with the dark matter sector 
and follows a common conservation law. In this work the authors have shown that the 
interacting EDE infact produce either eternal acceleration or deceleration, unless the 
dark energy density is modidifed with an additonal constant. 

6. On considering interaction between the cosmic components 
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